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Preface 


This  study  was  an  analytical  investigation  of  the  laminar,  steady, 
two-dimensional  extended  boundary  layer  equations  at  low  Reynolds 
number.  The  extended  boundary  layer  equations  are  the  Navier-Stokes 
equations  with  the  assumption  that  the  normal  pressure  gradient  is 
zero,  such  that  the  normal  component  of  the  momentum  equation  need  not 
be  solved.  This  differs  from  the  boundary  layer  equation  model  in  that 
the  streamwise  diffusion  term  is  not  neglected.  An  "extended  momentum 
integral  equation"  was  developed  from  the  continuity  equation  and  the 
streamwise  momentum  equation  and  applied  to  three  flow  problems  - flat 
plate,  stagnation  point,  circular  cylinder  - using  the  Karman- 
Pohlhausen  method. 

I would  especially  like  to  thank  my  thesis  advisory  committee. 
Professor  Richard  Merz  and  Professor  Harold  Wright,  for  their 
guidance  throughout  this  task.  I also  would  like  to  thank  my  faculty 
advisor.  Professor  Stephen  Koob,  for  his  guidance  and  much  needed  help 
in  this  study.  Finally,  I would  like  to  thank  my  wife,  Karen,  for  her 
tinderstanding  aind  encouragement  during  this  project. 

Robert  D.  Behr 
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Abstract 

An  analytical  study  of  the  application  of  the  Karman-Pohlhausen 
method  to  the  extended  boundary  layer  equations  at  low  Reynolds  number 
was  made.  The  extended  boixndary  layer  equations  were  the  incompressible 
Navier-Stokes  equations  with  the  assumption  of  zero  normal  pressure 
gradient.  A comparison  was  made  between  the  solutions  for  the  extended 
boundary  layer  equations  and  the  boundary  layer  equations  at  several 
Reynolds  numbers  for  flow  over  a flat  plate,  flow  near  a stagnation 
point,  and  flow  over  a circular  cylinder.  Favorable  results  were 
achieved  for  the  flat  plate  case  while  less  than  satisfactory  results 
were  achieved  for  the  circular  cylinder.  The  point  of  separation  was 
not  found.  Any  further  study  in  this  area  should  consider  using  ij 

cylindrical  coordinates  for  the  cylinder  solution  and  examining  the 
effect  of  the  pressure  gradient  normal  to  the  body. 
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! APPLICATION  OF  THE  KARMAN-POHLHAUSEN  METHOD 

TO  THE  EXTENDED  BOUNDARY  LAYER  EQUATIONS 

I.  Introduction 

Backgroimd 

In  the  derivation  of  the  boundary  layer  equations  from  the  Navier- 
Stokes  equations  it  is  assumed  that  the  pressure  /gradient  normal  to  the 
surface  is  zero  such  that  the  pressure  at  the  outer  edge  of  the  boundary 
layer  is  impressed  across  the  boundary  layer.  Furthermore,  the  velocity 
component  parallel  to  the  surface,  u,  becomes  equal  to  the  velocity  of 
the  outer  flow,  U^,  at  the  outer  edge  of  the  boimdar^'  layer.  The 
resulting  equations  are  known  as  Prandtl's  boimdaiy  layer  equations. 

For  low  Reynolds  number  flows  and  in  many  other  situations  (particularly 
large  adverse  pressure  gradients) , these  equations  do  not  adequately 
model  the  fluid  physics  and  the  full  Navier-Stokes  equations  should  be 
solved.  However,  due  to  their  coupled  non-linear  nature,  they  remain 
very  difficult  to  solve.  Therefore,  the  extended  boundary  layer 
equations  are  introduced  in  an  attempt  to  retain  some  of  the  physical 
modeling  accuracy  of  the  Navier-Stokes  equations  along  with  the 
simplicity  and  ease  of  solution  of  the  boundax-y  layer  equations.  The 
extended  boundary  layer  equations  are  the  two-dimensional  continuity 
and  streamwise  momentum  equations,  namely, 

= 0 (1) 

-f-VUy)  - -p^  tXX.  Uy^)  (2) 

irtiere  the  streamwise  viscous  diffusion  term,MU^,  is  not  neglected  as 
* in  the  boundary  layer  approximation.  This  work  addresses  an  "extended 
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momentum  Integral  equation"  for  the  extended  boundary  layer  equations 
; analogous  to  the  classical  momentum  integral  equation  for  the  boundary 

' layer  equations  and  solves  the  equation  for  several  external  pressure 

distributions  and  Rejrnolds  numbers. 

The  problem  of  flow  over  a flat  plate  has  previously  been  solved 
using  the  extended  boundary  layer  equations  by  Koob  and  Abbott  (Ref  2: 
64-7?) . This  previous  work  was  used  as  a starting  point  once  the 
momentum  integral  equation  was  developed  and  served  as  a procedural 
check  prior  to  examining  the  more  difficult  problems. 

Objectives 

The  objectives  of  this  work  were  as  follows: 

1.  Develop  a momentum  integral  equation  for  the  extended 
[ boundary  layer  equations. 

I 2.  Solve  the  extended  momentum  integral  eq^uation  numerically 

I for  several  Reynolds  numbers  and  the  following  external  pressure 

[ 

^ distributions : 

’ a.  Flow  over  a flat  plate. 

[ b.  Stagnation  point  flow. 

! c.  Flow  over  a circular  cylinder. 

[ 3*  Compare  the  extended  boundary  layer  and  Prandtl  boundary 

layer  results . 

Approach 

The  extended  momentiam  integral  equation  was  developed  with  the 
following  assumptions: 

1.  Flow  was  two-dimensional,  laminar,  incompressible,  and 

steady. 

2.  Pressure  was  only  a function  of  distance,  x,  along  the 

u 


1 


i] 


2 
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3t  Reynolds  number  was  small. 

4.  Interaction  between  the  external  flow  and  boundary  layer 


^ ; 

! 

i 

was  negligible . 

The  extended  momentum  integral  equation  was  then  non-dimensionalized 
! and  polynomials  of  various  degrees  were  assumed  for  the  velocity  profile 

to  transform  the  extended  momentum  integral  equation  from  a differential 
equation  in  four  dependent  variables  to  a second  order,  ordinary 
differential  equation  in  one  dependent  variable.  The  dependent 
variable  used  was  the  non-dimensional  displacement  thickness,  S, . A 
fourth  order  Runge-Kutta  computer  program  was  then  used  to  solve  the 
differential  equation  with  specified  initial  conditions.  Boundary 
layer  and  extended  boundary  layer  results  were  compared.  It  was 
desired  to  show  that  the  extended  boundary  layer  equations  would  be  a 
more  accurate  model  of  fluid  flows  than  the  boundary  layer  equations, 

I especially  in  adverse  pressure  gradients. 


II . The  Approximate  Solution  Method 


Development  of  the  Extended  Momentum  Integral  Equation 

Since  the  full  Navier-Stokes  equations  are  very  difficxilt  to  solve 
and  the  boundary  layer  equations  do  not  apply  near  and  beyond  separation, 
the  extended  boundary  layer  equations,  Eq  1 and  Eq  2,  were  proposed  as 
a new  model  to  meet  accuracy  and  ease  of  solution  requirements. 

Basically,  the  extended  boundary  layer  equations  were  the  two-dimensional 
Incompressible,  steady,  Navier-Stokes  equations  with  the  assumption  of  a 
zero  normal  pressure  gradient.  From  the  two-dimensional  continuity 
equation,  Eq  1,  and  the  streanwise  momentum  equation,  Eq  2,  an  "extended 
momentum  integral  equation"  was  developed.  At  the  edge  of  the  boundary 
layer,  where  the  boundary  layer  has  merged  with  the  external  flow,  Eq  2 
had  the  form: 

/I4l4^  = -Px  (3) 

where  is  the  velocity  at  the  edge  of  the  boundary  layer.  Substituting 
the  relationship  for  p^  in  Eq  3 into  Eq  2 resulted  in  the  following 
equation : 

/*  U u.^  + / iru.y  ^ file.  + M Uyy)  (4) 

The  same  procedure  described  in  Schlichting  (Ref  145)  for  the 
development  of  the  boundary  layer  momentum  integral  equation  was  then 
followed  tc  develop  the  extended  momentum  integral  equation. 

With  the  inclusion  of  the  streamwise  diffusion  term,^u^,  the 
final  equation  differed  from  that  presented  in  Schlichting.  Two 
additional  terms  appeared  from  integration  of  the  term,^u^.  The 
form  of  this  term  was  as  follows  (after  Eq  4 had  been  divided  through 


r ■■ 


"by  f ajid  the  term  Vu  introduced  to  aid  the  integration)! 

®xx 


(5) 


Thus,  the  extended  momentum  Integral  equation  was  developed  to  te 


where  the  shear  stress,  is  defined  by 

M,Uy\ 


XX 


(6) 


^ as  'y*0  ^ -^1  f I 

the  displacement  thickness,  S^t  is  defined  by 

* 

Mb.  S ~ J L^e.~  u)  dy 


(7) 


(8) 


and  the  momentum  thickness,  9,  is  defined  by 

^ H 

o = J uCUe.' 


(9) 


Since  Eq  6 was  in  a dimensional  fom,  it  was  advantageous  to 
non-dimensionalize  the  equation  with  respect  to  the  free  stream  velocity, 
Uj,  , and  a characteristic  length,  JL  . The  following  non-dimensional 
variables  were  used: 


V~Ue./Uco 

\'VJi 

Making  the  proper  substitutions  and  defining  Reynolds  number  as 


n . 

Ke  M. 


Eq  6 became 


= (Sx),  + (2Sx  ^-S.)  ^ + cfj  ^ 


'f 


(10) 


where  h in  Eq  6 has  been  replaced  by  the  boundary  layer  thickness,  S . 


I 


As  seen,  Eq  10  vjas  still  a second  order  differential  equation  with 
four  dependent  variables  - S , S,,  u^.  | y_Q.  The  Karman-Pohlhausen 

method  of  assuming  a polynomial  of  various  degrees  for  the  velocity 
profile  was  used  to  establish  three  additional  relationships  among 
these  dependent  variables. 

Assmned  Velocity  Profiles 

The  Karmain-Pohlhausen  method  has  been  demonstrated  to  be  a useful 
approximate  method  for  solving  the  boundary  layer  equations.  The  sane 
techniques  and  conditions  were  applied  to  the  extended  momentum  integral 
equation,  Eq  10,  assuming  polynomials  of  the  first  through  fourth  degree 
for  the  velocity  profile.  The  first  and  second  degree  polynomials  were 
used  primarily  to  establish  those  procedures  in  integration  which  would 
provide  the  best  results  for  the  other  velocity  profiles. 

When  choosing  the  velocity  profile,  consideration  was  given  to 
requiring  certain  properties  or  conditions  to  exist.  These  were  taken 
from  Schlichting  (Ref  3^  193)  S'Hd  were  as  follows; 

1.  No-slip  condition  at  the  wall. 

u(x,o)*o  (11) 

2.  Continuity  at  the  point  where  the  boundary  layer  meets 
the  potential  flow. 

u.(x,S)*Ue  (12) 

3.  Zero  shear  stress  at  the  boundary  layer  edge, 

Uy(x,S)^0  (13) 

4.  Satisfy  the  streamwlse  momentum  equation  at  the  body 

surface . i 

l/u.,y(x,o)  ~Uc  (14) 
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5.  Satisfy  the  inviscid  (potential)  equation  at  the  boundary 
layer  edge. 

(15) 

With  Eq  11  through  Eq  15  serving  as  the  boundary  conditions,  a 
velocity  profile  was  assumed  in  the  form  of  a polynomial  in  terms  of 
the  non-dimensional  distance  from  the  surface 

1 • 

A general  expression  for  this  polynomial  was  as  follows: 

* a •»-  b)j  . ..  (16) 

where  a,b,c,d,e, . . . were  functions  of  x to  be  evaluated  from  the 
boundary  conditions.  The  boundary  conditions  used  depended  on  the 
degree  of  the  polynomial  selected  and  the  judgement  of  the  author. 

Table  I gives  a summary  of  the  velocity  profiles  examined  in  this 
study  with  the  boundary  conditions  used  and  three  of  the  dependent 
variables  defined  as  functions  of  the  boundaiy  layer  thickness,  S , 
and  characteristic  length,  Jt.  The  dimensionless  quantities  J amd  K 
appeared  as  a result  of  applying  the  boundary  conditions  of  Eq  14  and 
Eq  15.  In  equation  form. 


J - * ({4/  Re  Vf 

(17) 

K - - cw)*^ 

(18) 

Differential  Equation  in  One  Dependent  Variable.  S. 

With  the  velocity  profile  assumed,  the  four  dependent  variables  - 
S*  S,»  I y=o  ” were  related.  Therefore  it  was 

possible  to  represent  three  of  the  dependent  variables  in  terms  of 
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of  Velocity  Profile  Properties  and  Results 


the  fourth  thereby  reducing  the  number  of  dependent  variables  in  the 
equation  to  one.  The  variable  selected  for  use  vfas  the  non-dimensional 
displacement  thickness , S, , Making  the  proper  substitutions  into  Eq  10 
for  5^  and  u^j^^  for  each  respective  velocity  profile  resulted  in  the 
second  order  differential  equations  listed  in  Table  I as  Eq  19  to  Eq  22. 
The  symbol  ' infers  differentiation  with  respect  to  ^ . 

Numerical  Integration 

The  solution  to  the  second  order  differential  equation  was  obtained 
by  nximerical  integration  using  a fourth  order  Runge-Kutta  computer 
program  as  discussed  in  Ketter  and  Prawel  (Ref  1;  268-277) • The  step 
sizes  used  in  the  program  were  adjusted  according  to  Table  II.  Any 
additional  constraints  or  requirements  on  the  integration  will  be 
discussed  in  the  appropriate  sections. 


Table  II 


Integration  Step  Size  and  Corresponding  Range  of  ^ 


Step  Size 

Range 

* ?/io 

0.000  ^ C 0.01 

0.005 

0.010  ^ ^ < 1.0 

0.010 

1.000  ^ < 5.0 

0.050 

5.0  ^ < 10.0 

0.100 

10.0  ^^<100.0 

0.500 

100.0  ^^<  1000.0 

9 


Ill . Flat  Plate  Solution 


The  flat  plate  problem  tdiich  has  no  external  pressxire  gradient 
was  solved  by  Koob  and  Abbott  (Ref  2),  In  their  study,  the  extended 
boundary  layer  equations  were  solved  numerically  by  requiring  the 
numerical  solution  to  form  a bridge  between  the  Stokes  solution, 
which  is  valid  near  the  leading  edge,  and  the  Prandtl  approximation 
solution,  which  is  valid  for  large  Reynolds  numbers.  An  Initial  S,  was 
chosen  sufficiently  close  to  the  leading  edge  so  that  the  Stokes 
equation  and  extended  boundary  layer  equation  were  equivalent.  Using 
the  Stokes  approximation  solution  to  determine  starting  conditions, 
the  extended  boundary  layer  differential  equation  was  numerically 
integrated  in  an  attempt  to  reach  a point  where  the  Blasius  solution 
would  be  realized.  However,  the  integration  could  not  be  carried  far 
enough  for  the  Blasius  character  to  be  apparent  and  an  asymptotic 
exp2insion  was  used  for  ^ greater  than  45  to  compute  S,  and  eventually 
reach  the  Blasius  solution.  The  results  from  integrating  the  extended 
boundary  layer  equations  were  in  agreement  with  the  Stokes  solution  for 
^ less  than  0.25  ajid  with  the  Prandtl  solution  for  ^ greater  than  100. 
These  results  were  used  to  evaluate  the  procedures  applied  in  this  study 
and  to  provide  a comparison  of  the  extended  boundary  layer  solutions  for 
this  simple  case  before  undertaking  the  more  difficult  problems  of  the 
stagnation  point  and  circular  cylinder.  Thus,  the  procedures  to  be 
followed  here  were  analogous  to  those  in  Koob  and  Abbott  (Ref  2). 

The  applicable  differential  equations  came  from  Eq  19  through 
Eq  22  with  substitution  made  for  the  non-dimensional  external  pressure 
distribution  and  its  derivatives.  For  the  flat  plate  problem,  U 
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equaled  one  and  all  its  derivatives  were  zero.  Thus,  the  variables  J 
and  K,  defined  by  Eq  17  and  Eq  18,  respectively,  were  also  zero.  The 
equations  were  significantly  simplified  and  the  techniques  employed  by 
Koob  and  Abbott  (Ref  2:  66-70)  were  applied  for  all  four  velocity 
profiles.  To  demonstrate  the  techniques,  as  they  apply  to  the 
differential  equations,  the  third  degree  velocity  profile  equation 
will  be  developed  here. 

Beginning  with  Eq  21  the  following  substitutions  were  made; 

V-l 

V'=U*=6) 

J - J' » o 

Ibe  result  was 


d,  ® d,  1 


(23) 


The  corresponding  boundary  layer  equation  was  derived  by  neglecting 
the  last  two  terms  in  Eq  10  or  by  setting  equal  to  zero  in  Eq  23 

g,  - (24) 

Since  the  boundary  layer  equations  were  usually  transformed  to  be 
independent  of  Reynolds  number,  a transformation  was  performed  to 
remove  the  Reynolds  number  from  Eq  23  and  Eq  24  by  assuming 

X “ s. 

Equation  23  became 

S,  Re  » (ssjRe  ~ 7tf,  ^ 

Equating  the  exponents  of  Reynolds  number  gave 


11 


from  which  the  values  of  f and  ^ were  determined  to  be  unity 


\ 


I 

1 

I 


This  transformation  resulted  in  Eq  23  and  Eq  24  being  Reynolds  number 
independent.  The  procedure  was  repeated  for  the  other  profiles  and 
yielded  the  same  results  for  f and  ^ . Table  III  summarizes  the 
extended  botandary  layer  and  boundary  layer  equations  for  the  various 
velocity  profiles.  The  bars  on  S,  have  been  dropped  for  convenience. 

Table  III 

Equations  for  Flow  Over  a Flat  Plate 


Velocity  Profile 

Extended  Boundary 
Layer  Ea nation 

Boundary 
Layer  Equation 

Stokes  Equation 

First  degree 

s:-s.u)-^ 

i:'= 

a.  zi. 

s - J- 
0| 

Second  degree 

If 

5,  - 

Third  degree 

s!" ' - ih, 

1 3/S^ 

S,  - zoBSi 

Fourth  degree 

5."*  5, 

c'  = .gfeZ, 

0,  3 70^, 

c" 2- 

di  " s-b, 

The  extended  bovmdary  layer  equation  was  singular  for  the  first 
of  two  boundary  conditions 

O aT  ^ -o  (25) 

S,  — ► BIasius  Limit  at  ^ ► 00  (26) 

However,  two  solutions  which  were  valid  near  ^ equal  to  zero  and  ^ 
approaching  infinity  were  determined  from  the  Stokes  equation  and  the 
boundary  layer  equation  in  Table  III . The  extended  boundarj-  layer 
equation  was  solved  numerically  to  form  a bridge  between  these  two 
solutions . 

Again  using  the  third  degree  velocity  profile  as  an  example,  the 

12 


I 


I 


Stokes  equation  from  Table  III  was 

■=  " 71?,  (27) 

Equation  2?  was  analytically  integrated  by  introducing  the  variable 


such  that , 

p = ^C--|  Ih  8,  (28) 

Integration  of  Eq  28  yielded 

^ e^^erfc  C-lnS,  ) (29) 

vhere  C was  a constant  of  integration  which  must  be  evaluated. 

Following  the  procedures  of  Koob  and  Abbott  (Ref  2;  71),  Eq  28  and 
Eq  29  were  used  to  determine  the  initial  conditions  for  the  integration 
of  the  extended  boundary  layer  equation,  Eq  23. 

1.  Pick  values  for  the  integration  constant,  C,  and  initial  i',. 

/ - 

2.  Solve  Eq  28  and  Eq  29  for  5,  and  y , respectively. 

3.  Integrate  Eq  23  with  these  initial  conditions, 

4.  Adjust  value  of  C until  the  extended  boundary  layer 
solution  approaches  the  boundary  layer  solution  at  ^ 
approaching  infinity. 

-12 

The  Initial  value  of  S,  selected  was  10  . During  the  Integration 

process  if  S,  became  less  than  zero  or  6,  became  greater  than  zero,  the 
Integration  was  stopped  since  this  would  be  in  violation  of  the  asymptotic 
behavior  of  S,  (See  Figure  1),  A new  value  of  C was  then  determined  and 
the  procedure  repeated.  It  was  during  this  procedure  that  the  range 
and  size  of  the  integration  step  were  refined  (See  Table  II).  The  value 


of  C differed  for  each  velocity  profile?  however,  the  limits  in  the 


Table  V depicts  the  results  of  the  extended  boundary  layer  and 
boundary  layer  approximations  for  each  case  and  the  per  cent  difference 
as  compared  to  each  other.  This  per  cent  difference  was  computed  by 
taking  the  difference  between  the  values  of  the  two  solutions  and 
dividing  by  the  extended  boundary  layer  solution  value.  As  shown,  the 
extended  boundary  layer  and  boundary  layer  solutions  agreed  within  five 
I)er  cent  of  one  another  at  ^ =50,  but  for  smaller  ^ the  error  increased 
rapidly.  These  same  results  were  obtained  by  Koob  and  Abbott  (Ref  2:  75) 
using  the  Blasius  profile. 


Table  V 


I 

I 

Results  and  Comparison  of  the  Extended  Boundary/  Layer 
and  Boundary  layer  Solutions  for  a Flat  Plate 


First  Degree  Velocity  Profile 


Si  Extended 

s, 

6 

Boundary  Layer 

Boundary  Layer 

% Difference 

0.001 

0.002588 

0.05477 

2016. 

0.01 

0.02138 

0.1732 

710. 

0.1 

0.1606 

0.5^77 

241. 

1.0 

1.003 

1.732 

72.7 

10.0 

4.667 

5.477 

17.35 

50.0 

11.620 

12.247 

5.39 

Second  Degree 

Velocity  Profile 

8,  Extended 

s, 

5 

Boundaiy  Layer 

0.001 

0.002950 

0.05774 

1857. 

0.01 

0.02424 

0.1826 

653. 

0.1 

0.1801 

0.5774 

220. 

1.0 

1.102 

1.826 

65.7 

10.0 

5.007 

5.774 

15.30 

50.0 

12.330 

12.910 

4.70 

Third  Degree 

Velocity  Profile 

8,  Extended 

6. 

■■■■ 

t _ 

Boundary  Layer 

Boundary  Layer 

0.001 

0.002724 

0.05503 

1920. 

0.01 

0.02243 

0.1740 

675. 

0.1 

0.1675 

0.5503 

229. 

1.0 

1.033 

1.740 

68.4 

10.0 

4.740 

5.503 

16.10 

50.0 

11.723 

12.306 

4.97 

Fourth  Degree  Velocity  Profile 


5 

S,  Extended 
Boundary  Layer 

6. 

Boundary  Layer 

0.001 

0.002803 

0.05536 

1875. 

0.01 

0.02304 

0.1751 

660. 

0.1 

0.1715 

0.5536 

222. 

1.0 

1.052 

1.751 

66.4 

10.0 

4.792 

5.536 

15.52 

50.0 

11.816 

12.379 

4.76 
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IV.  stagnation  Point  Solution 

Introduction 

The  boundary  layer  flow  around  a shaped  body  generally  proceeds 
from  a stagnation  point  to  a point  of  separation  and  beyond.  Therefore, 
the  staignatlon  point  problem  was  considered  in  order  to  establish  the 
proper  initial  conditions  for  the  circular  cylinder  problem  which  is  seen 
as  the  principal  test  of  the  new  method.  The  discussion  was  limited  to 
the  assumed  third  and  fourth  degree  polynomial  velocity  profiles  since 
the  lower  degree  profile  shapes  were  independent  of  pressure  gradient. 
Exact  Solution 

An  exact  solution  to  the  complete  Kavier-Stokes  equations  was 
fovind  by  K.  Hiemenz  (Ref  3^  87-90).  The  potential  flow  was  defined 
by 

ax 
Vt  = -a.y 

► where  "a"  denotes  a constant.  As  shown  in  Figure  2,  the  potential 

[ 

[ flow,  perpendicular  to  a flat  wall,  impacts  and  leaves  the  area  of 

I impact  in  two  directions  along  the  x-axis. 
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For  the  potential  flow,  the  pressure  was  given  by  Bernoulli's  equation 


Po  - p * i /*  ( 14^+  Ve  j 

A similarity  solution  was  assumed 

U.=  Xf\y) 

V*  --fCy) 

P.-P  = ( X^l-Fcy)) 

which  identically  satisfied  the  continuity  equation.  The  differential 
equations  for  f and  F from  the  Navier-Stokes  equations  were 

(30) 

■ff'  = -i:a  F'-Vf"  (31) 

The  boundary  conditions  at  the  wall,  where  u=v=0,  and  at  a large 

distance  from  the  wall,  where  u=U  =ax,  determined  the  boundary  conditions 

e 

on  f and  F . 

f F » 0 of  y~0  (32) 

f CL  odr  CO  (33) 

A solution  for  f was  determined  first  from  Eq  30  by  solving  it 
numerically  since  it  could  not  be  solved  in  closed  form.  Equation  30 
was  non-dimensionalized  by  assuming 

Vj- 

where 

o<  - 

This  led  to  the  differential  equation 

»***  . ^ -'A 

<p  i-  <p^  - <p 
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where  ' denotes  differentiation  with  respect  to  . The  boundary 
conditions  were  converted  to 

<p'  - O at  yj  ~ O 

<p‘  ^ t At  yj  >=  oa 

Thus,  the  velocity  component,  u,  became 

fre"  (3^) 

As  increased,  increased  linearly  at  first  and  then  asymptotically 
approached  one.  At  = 2.4  the  value  of  was  approximately  0.99 
or  within  one  per  cent  of  the  final  value.  If  the  boundary  layer  edge 

is  defined  to  be  the  value  of  y = S v^ere  u = 0.99  U^,  then 

S “ (35) 

where  was  that  value  of  where  u » 0.99  U^*  A complete  table  of 
y\  t ^ and  (p'  is  found  in  Schlichting  (Ref  3*  90) . 

Approximate  Solution  from  the  Extended  Momentum  Integral  Equation 

The  flow  in  the  vicinity  of  the  staignation  point  was  characterized 
by  the  external  velocity  distribution 

zVX 

or  in  dimensionless  coordinates 

U = (36) 

With  Eq  36  applied  to  the  appropriate  differential  equations  in  Table  I , 
the  stagnation  point  equations  were  developed  and  recorded  in  Table  VI. 
Removal  of  Reynolds  number  was  not  possible  since  Reynolds  number  is 
contained  in  the  definition  of  J.  To  numerically  Integrate,  it  was 
necessary  to  determine  initial  values  for  S,,  8, , £ 1 and  the  quantities 
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J and  K.  For  the  exact  solution,  Eq  35»  S did  not  vary  with  x or  ^ ; 
therefore,  S,  was  also  constant  and  one  of  the  initial  conditions 
prescribed  was  that 

s/  = 0 at  ^ - 0.0001 

irtiere  the  value  of  ^ was  chosen  arbitrarily  small.  The  value  of  K 
was  zero  since  from  Eq  36  is  zero.  From  Eq  I7,  the  J equation  for 


stagnation  point  flow  became 


(41) 


Using  the  expressions  in  Table  I for  the  displacement  thickness,  , 
a relationship  between  J and  was  derived  for  each  profile.  These 
relationships  were 

X - ^^8  S* 

vi  “ 


(42) 


for  the  third  degree  velocity  profile  and 

T - ^00 

0 “ 


(43) 


for  the  fourth  degree  velocity  profile.  Thus,  for  a given  Reynolds 
number,  the  values  of  J and  S,  must  satisfy  Eq  42  or  Eq  43.  The 
possibility  existed  that  J would  have  three  values  (roots)  for  a given 
S,  and  Reynolds  number. 

Determination  of  Initial  Values  for  J and  6. 

With  S,  constant  for  the  stagnation  point  flow,  the  following 
condition  was  implied: 

at  (44) 

From  Eq  41  or  the  complementary  equations,  Eq  42  or  Eq  43,  the  value  of 
J was  also  constant  for  a given  Reynolds  number.  Thus, 

j-'«  o (45) 
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Using  the  conditions  of  Eq  ^ and  Eq  45,  the  extended  boundary  layer 
equations  in  Table  VI  were  reduced  to  functions  of  S,  , J and  Reynolds 
niimber.  These  new  equations  were  then  used  with  Eq  42  or  Eq  43  to 
solve  for  J and  then  S,  . For  the  fourth  degree  velocity  profile  Eq  44 
and  Eq  45  were  imposed  on  Eq  39  with  the  following  results 

MjJ- 

or,  after  rearranging, 

ZO  $,  Re  = 33^,  - 71^^^  - ^ (2^) 

Substituting  Eq  46  into  Eq  43  gave 

J “ (^7) 

which  is  a third  degree  equation  for  J.  The  roots  of  J for  both  the 
fourth  degree  and  third  degree  velocity  proflZes  are  given  in  Table  "/II . 

Table  VII 

Values  of  J and  ifRe  Satisfying;  Equation  44 
and  Equation  42  or  Equation  43 


Third  Degree  Velc 

J 

)city  Profile  (42) 
S*Re 

Fourth  Degree  Ve] 

J 

Locity  Profile  (43) 
SfRe 

5.8653997  + 
2.1656705  i 

85.814115  + 

0.0126237  i 

7.0523231 

0.20519518 

5.8653997  - 
2.1656705  i 

85.814115  - 

0.0126237  1 

I7.8O3257 

0.2046889 

-32.230799 

-17.648133 

-72.255580 

-29.402177 

From  Eq  41,  J must  be  real  and  positive;  therefore,  a solution  for  the 
third  degree  velocity  profile  that  meets  the  required  conditions  for 
the  stagnation  point  was  not  possible.  The  reason  for  this  is  not 
understood  at  this  time.  For  the  fourth  degree  velocity  profile  two 
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I 

j 


roots  for  J seemed  plausible: 


J,  - Z0S2323I 
J,  = /7.S03257 

From  the  exact  solution  discussed  previously  and  in  Schlichting 
(Ref  3*  89)  a relationship  was  developed  between  J and  Equation 

41  was  rewritten  as  follows; 

J - 2.  (S/i)  = 2 (S/l)  = S (■^) 

:diere  "a”  was  transformed  from 

a-  Oe/x 

to 

0.*  2 

Applying  Eq  35  gave 

J*  (48) 

Corresponding  values  of  were  computed  for  and  Jg  from  Eq  48.  The 
»j  values  in  Table  5*1  of  Schlichting  (Ref  3s  90)  were  redefined  to 
correspond  to  »]  = y/S  used  in  the  assumed  velocity  profiles  by 
dividing  the  rj  in  the  table  by  the  computed  :jj  for  the  respective  J 
value.  With  the  exact  solution  now  defined  in  the  same  terms  as  the 
approximate  solution,  a comparison  was  made  between  the  two  solutions 
based  on  the  J value.  These  comparisons  are  presented  in  Figures  3 
and  4. 

Both  and  were  tried  as  the  J value  for  the  stagnation  point. 
For  a given  Reynolds  number,  S,  was  derived  from  Eq  43  (See  Table  VIl). 
Thus,  the  required  Initial  conditions  were 


i 

■ 


I 

« 
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■ 
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S>0 

J,-  7.0S23231 

s,  ='^0.20511518/^ 

and  the  integration  of  Eq  39  was  hegm. 

Results 

Throughout  the  range  of  integration,  S,  remained  constant  or 
nearly  so  for  both  the  extended  botindary  layer  and  boimdary  layer 
solutions.  This  demonstrated  that  the  equations  satisfied  the  conditions 
at  or  near  xhe  stagnation  point.  With  the  starting  conditions  defined, 
the  stagnation  point  results  were  used  as  initial  conditions  for  the 
circular  cylinder  flow. 


§.'-0 

J,-  17.803257 
S,  -'^0.2oH(,BByiSe. 


i 

1 


\ t 
i ■ 

I i i 

.J.  * 


I 


il 


1. 


and  Fourth  Degree  Profile  Solution  with  J = 7.052323I 
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Figure  4.  Comparison  of  Velocity  Distributions  for  Exact  Solution 
and  Fourth  Degree  Profile  Solution  with  J = 17.803257 
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V.  Circular  Cylinder  Solution 


Introduction 

The  flow  around  a circular  cylinder  Is  a classical  case  for 
testing  viscous  equation  solution  methods  as  It  permits  the  study 
of  a boundary  layer  from  a steignatlon  point  throiigh  a separation  point. 
The  flow  Is  characterized  by  the  external  velocity  distribution 

U = 2 $in(X/R) 


or  In  dimensionless  coordinates, 


U = 2 sin  ^ 


where  ^ Is  the  angle  In  radians  measured  from  the  stagnation  point 
and  R Is  the  cylinder  radius  now  used  as  the  characteristic  length. 
The  discussion  is  limited  to  the  solution  for  the  assumed  fourth 
degree  velocity  profile  since  the  conditions  at  the  stagnation 
point  for  the  third  degree  velocity  profile  were  undetermined.  The 
equations  for  the  circular  cylinder  are  listed  in  Table  VIII. 

With  the  given  pressure  distribution  the  values  of  J and  K were 
as  follows: 


T siooS,  Rc-  ^ 

^ - /oo  S, — ^ 

^ O-Vizi-'y^y  2 Re.  cos  ^ 

The  initial  conditions  imposed  for  this  case  were  taken  from  the 


(^9) 

(50) 


stagnation  point  results 


J - Z 052323  I or 


at  ^ 0.000\ 

/Z  S032  5 7 


where  the  initial  value  of  ^ was  arbitrarily  selected.  The  initial 


Circular  Cvllnder  Equations  for  Fourth  Deirree  Velocity  Profile 


values  of  K and  S,  were  then  computed  based  on  the  selected  J value 
£uid  Reynolds  number. 

Results 

With  the  initial  values  computed,  the  integration  was  begun  and 

4 

proceeded  until  S,  became  greater  than  10  or  until  ^ equaled  J.lk, 
indicating  the  aft  portion  of  the  cylinder  had  been  reached.  The 
extended  boundary  layer  and  boundary  layer  solutions  were  integrated  at 
Reynolds  numbers  of  1,  10,  100,  and  1000.  For  a given  Rej-nolds  number, 
the  extended  boundary  layer  and  boundary  layer  solutions  began  from 
nearly  the  same  value  of  S,  . However,  b,  was  not  zero  or  nearly  zero 
for  the  extended  boundary  layer  solution  near  the  stagnation  point  as 
was  the  case  for  the  stagnation  point  flow.  This  was  the  result  of 
the  term,  U"/U,  in  Eq  22  no  longer  being  zero.  This  term  caused  S, 
to  have  a large  negative  value  and  the  extended  boundary  layer  solution 
(€,)  decayed  rather  quickly  for  the  fourth  degree  velocity  profile  as 
depicted  in  Figure  5*  Meanwhile,  the  boundary  layer  solution  had  S, 
increasing  gradually  until  a certain  angle  from  the  stagnation  point, 

(f>  , was  approached,  at  which  time  the  S,  increased  rapidly  and  exceeded 

4 . 

10  . The  higher  the  Reynolds  number,  the  higher  was  the  value  of  (p 

where  6,  increased  rapidly.  The  point  of  rapid  S,  increase  was  assumed 

to  be  where  separation  occurs.  As  shown  in  Figure  5»  "the  angle  ^ was 

in  the  ramge  of  I30  degrees  to  I50  degrees  for  the  Reynolds  numbers 

tested  rather  than  approximately  82  degrees  as  previous  studies 

have  shown  (Ref  ’}>%  202).  Thus,  the  present  results  are  not  understood. 

A possible  reason  for  this  conflict  could  be  the  coordinate  system 
used.  A rectangular  coordinate  system  does  not  account  for  the  curvature 
of  the  cylinder.  This  could  Impact  the  solution  for  the  extended 


I 

J 
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bovindary  layer  sind  boundary  layer  approaches  dramatically.  Perhaps 
it  would  have  been  more  appropriate  to  have  used  a cylindrical 
coordinate  system  thereby  accounting  for  the  curvature. 


I 
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Fl^re  5»  Comparison  of  Extended  Boundary  Layer  and  Boundary  Layer 
Solutions  for  Reynolds  Numbers  of  1,  10,  and  100 


VI . Conclusions  and  Recommendations 


Conclusions 

The  conclusions  drawn  from  this  study  of  the  extended  boundary 
layer  equations  are  presented  as  follows: 

1.  A momentum  integral  equation  can  be  developed  which 
includes  the  streamwise  diffusion  term.  It  is  hypothesized  that  this 
is  a more  accurate  approach  than  the  boundary  layer  equations  and  a 
necessary  approach  for  boimdary  layer  flows  near  separation. 

2.  The  results  obtained  for  the  flow  over  a flat  plate  were 
in  agreement  with  those  presented  in  other  studies, 

3.  The  results  obtained  for  the  flow  at  the  stagnation 
point  and  over  a circular  cylinder  were  less  than  satisfactory.  The 
extended  boundary  layer  and  boundary  layer  solutions  did  not  agree 
for  the  flow  over  a cylinder  at  any  Reynolds  number. 

k.  The  method  applied  did  not  give  a clear  indication  as 

to  where  separation  occurred  on  the  cylinder  as  compared  to  previously 
demonstrated  methods. 

Recommendations 

The  following  recommencations  are  made  for  future  studies  in 
this  area: 

l,  The  coordinate  system  used  in  the  extended  boundary  layer 
equations  should  be  changed  to  a cylindrical  coordinate  system, 
appropriate  equations  developed,  and  the  same  tests  run  to  determine 
idiether  correct  results  could  be  obtained. 

2.  The  y-momentum  equation  should  be  solved  for  the  normal 
pressure  gradient  term  to  determine  whether  it  is  negligibly  small  as 
assumed. 
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An  analytical  study  of  the  application  of  the  Karman-Pohlhausen  method  to 
the  extended  boundary  layer  equations  at  low  Reynolds  number  was  made.  The 
extended  boundary  layer  equations  were  the  Incompressible  Navier-Stokes 
equations  with  the  assumption  of  zero  normal  pressure  gradient.  A comparison 
was  made  between  the  solutions  for  the  extended  boundary  layer  equations  and 
the  boundary  layer  equations  at  several  Reynolds  numbers  for  flow  over  a flat 
plate,  flow  near  a stagnation  point,  and  flow  over  a circular  cylinder. 
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Favorable  results  were  achieved  for  the  flat  plate  case  while  less  than 
satisfactory  results  were  achieved  for  the  circular  cylinder.  The  point 
of  separation  was  not  found.  Any  further  study  in  this  area  should 
consider  using  cylindrical  coordinates  for  the  cylinder  solution  amd 
examining  the  effect  of  the  pressure  gradient  normal  to  the  body. 
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